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ABSTRACT 

We have obtained early-time photometry and spectroscopy of supernova (SN) 2013df in NGC 4414. The SN is 
clearly of Type lib, with notable similarities to SN 1993J. From its luminosity at secondary maximum light, it 
appears that less 56 Ni (<0.06 M 0 ) was synthesized in the SN 2013df explosion than was the case for the SNe lib 
1993J, 2008ax, and 201 ldh. Based on a comparison of the light curves, the SN 2013df progenitor must have been 
more extended in radius prior to explosion than the progenitor of SN 1993J. The total extinction for SN 2013df 
is estimated to be Ay — 0.30 mag. The metallicity at the SN location is likely to be solar. We have conducted 
Hubble Space Telescope ( HST) Target of Opportunity observations of the SN with the Wide Field Camera 3, and 
from a precise comparison of these new observations to archival HST observations of the host galaxy obtained 
14 yr prior to explosion, we have identified the progenitor of SN 2013df to be a yellow supergiant, somewhat hotter 
than a red supergiant progenitor for a normal Type 11-Plateau SN. From its observed spectral energy distribution, 
assuming that the light is dominated by one star, the progenitor had effective temperature T e g = 4250 ± 100 K and 
a bolometric luminosity Lb 0 i = 10 4,94±0 06 L G . This leads to an effective radius R e g = 545 ± 65 Rq. The star likely 
had an initial mass in the range of 13-17 M 0 ; however, if it was a member of an interacting binary system, detailed 
modeling of the system is required to estimate this mass more accurately. The progenitor star of SN 2013df appears 
to have been relatively similar to the progenitor of SN 1993J. 
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1. INTRODUCTION 

A core-collapse supernova (SN) is the final stage of evolution 
for stars with initial mass M mi > 8 M Q (e.g., Woosley et al. 
2002). The compact remnant of this explosion is thought to 
be either a neutron star or a black hole. Stars that reach the 
terminus of their evolution with most of their hydrogen-rich 
envelope intact — single red supergiants (RSGs) — produce Type 
II (specifically Type II-Plateau) SNe. As the envelope is stripped 
away, either through vigorous mass loss or via mass transfer in 
a binary system, the results are the hydrogen-free, yet helium- 
rich, Type lb SNe and the hydrogen-free and helium-poor (or 
helium-free) Type Ic SNe; see Filippenko (1997) for a review of 
SN classification. Intermediate to the SNe II and SNe lb are the 
Type lib, which retain a low-mass (<1 M G ) hydrogen envelope 
prior to explosion. SNe lib are intrinsically rare (~10%-11% 
of all core-collapse SNe; Smith et al. 2011; Li et al. 2011). The 
first identified case was SN 1987K (Filippenko 1988); well- 
studied examples are the nearby SN 1993J in Messier 81 (e.g., 
Richmond et al. 1996; Matheson et al. 2000) and SN 201 ldh 
in Messier 51 (e.g., Arcavi et al. 2011; Ergon et al. 2013). The 
progenitor channel for SN 1993J that has garnered most favor 
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over the years has been a massive interacting binary system 
(e.g., Podsiadlowski et al. 1993; Woosley et al. 1994; Maund 
et al. 2004; Stancliffe & Eldridge 2009). A similar model has 
been introduced for SN 201 ldh (Benvenuto et al. 2013). The 
progenitor is the mass donor via Roche lobe overflow. 

Spectroscopically, SNe lib exhibit at early times emission 
and absorption features due to H, resembling SNe II, which 
subsequently give way to the emergence of the He I absorption 
lines typical of SNe lb (e.g., Filippenko et al. 1993; Chomock 
et al. 2011), with broad H emission reappearing in the nebular 
phase (e.g., Filippenko et al. 1994; Matheson et al. 2000; 
Taubenberger et al. 2011; Shivvers et al. 2013). Although rarely 
seen owing to its very short duration early in the evolution of 
an SN, a rapid decline after an initial peak has been observed 
among SNe lib (e.g., Richmond et al. 1994, 1996; Roming et al. 
2009; Arcavi et al. 2011). This is interpreted as adiabatic cooling 
after the SN shock has broken out through the star’s surface; the 
duration of this cooling is governed primarily by the radius of 
the progenitor (Chevalier & Fransson 2008; Nakar & Sari 2010; 
Rabinak & Waxman 2011; Bersten et al. 2012). The SN light 
curves reach a secondary maximum, as a result of thermalization 
of the y-rays and positrons emitted during the radioactive decay 
of 56 Ni and 56 Co, followed by a smooth decline. The light-curve 
shapes resemble those of SNe lb (Arcavi et al. 2012). 
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Figure 1. Grayscale composite of rizJH images obtained with RATIR on 2013 
June 24, showing SN 2013df in NGC 4414. The position of the SN is indicated 
by the tick marks. North is up, east is to the left. 


As the number of all core-collapse SNe with directly identi- 
fied progenitor stars is very small (~20 at the time of this writ- 
ing), we have been extraordinarily fortunate that, up to this point, 
the progenitors of three SNe lib have been identified, including 
SN 1993J (Aldering et al. 1994; Van Dyk et al. 2002; Maund 
et al. 2004; Maund & Smartt 2009), SN 2008ax in NGC 4490 
(Crockett et al. 2008), and SN 201 ldh in Messier 51 (Maund 
et al. 2011; Van Dyk et al. 201 1). Each of these stars exhibits 
distinctly different properties, although all show indications of 
envelope stripping prior to explosion: the SN 1993J progeni- 
tor has been characterized as a K-type supergiant, with initial 
mass ~ 13-22 M 0 (Van Dyk et al. 2002; Maund et al. 2004); 
for SN 2008ax it was difficult to fit a single supergiant to the 
observed colors (Lietal. 2008; Crockett etal. 2008), and the pos- 
sible initial mass range for the progenitor is large, ~ 10-28 M Q ; 
and for SN 20 1 1 dh we now know that the ~6000 K yellow 
supergiant identified by Maund et al. (201 1 ) and Van Dyk et al. 
(2011) was indeed the star that exploded (Ergon et al. 2013; Van 
Dyk et al. 2013b), with initial mass ~12-15 M Q (Maund et al. 
2011; Murphy et al. 2011; Bersten et al. 2012). In addition, 
Ryder et al. (2006) may have detected at very late times the 
massive companion to the SN lib 2001ig in NGC 7424 (Silver- 
man et al. 2009). Furthermore, Chevalier & Soderberg (2010) 
have divided SN lib progenitors into those that are compact 
(radius ~10 n cm; e.g., SN 2008ax) and those that are extended 
(~10 13 cm; SN 1993J), although SN 201 ldh appears to be an 
intermediate case (Horesh et al. 2013). Each new example pro- 
vides us with an increased understanding of this SN subtype and 
of the massive stars that give rise to these explosions. 

In this paper we consider SN 2013df in NGC 4414, shown 
in Figure 1. It was discovered by Ciabattari et al. (2013) on 
June 7.87 and 8.83 (UT dates are used throughout). Cenko 
et al. (2013) provided spectroscopic confirmation on June 10.8 
that it is a Type II SN. The resemblance to SN 1993J at early 
times suggested to Cenko et al. that SN 2013df would evolve 
to be an SN lib. Here we present early-time photometric and 


spectroscopic observations of SN 2013df, which demonstrate 
that it is indeed an SN lib. In a telegram, we presented a 
preliminary identification of three progenitor candidates for the 
SN (Van Dyk et al. 2013a). We provide here a much better 
identification, through high-resolution imaging of the SN, and 
show that none of the three sources turned out actually to be at 
the SN position. We will characterize the nature of the probable 
progenitor. The host galaxy, NGC 4414, is a nearby, isolated, 
flocculent spiral galaxy, with an inclination of 55° (Vallejo 
et al. 2002). Thornley & Mundy (1997) measured a global star 
formation rate of 1.3 M Q yr _l , comparable to that of other Sc 
galaxies. SN 2013df is at a nuclear offset of 32" E, 14" N, along 
an outer spiral arm. The galaxy was also host to the likely Type 
la SN 1974G (e.g., Schaefer 1998). We adopt the Cepheid-based 
distance modulus /xo = 31.10 ± 0.05 mag (distance 16.6 ± 0.4 
Mpc) established by Freedman et al. (2001). 

2. OBSERVATIONS 
2.1. Early -time Photometry 

We have observed the SN using the 0.76 m Katzman Au- 
tomatic Imaging Telescope (KAIT; Filippenko et al. 2001) at 
Lick Observatory in BVRI between 2013 June 13.7 and July 
18.7. Point-spread function photometry was applied using the 
DAOPHOT (Stetson 1987) package from the IDL Astronomy 
User’s Library 1 1 . The instrumental magnitudes and colors of the 
SN were transformed to standard Johnson-Cousins BVRI using 
two stars in the SN field which are in the Sloan Digital Sky 
Survey (SDSS) catalog, and by following this prescription 12 to 
convert from the SDSS to the Johnson-Cousins system. 

Data in rizJH were also obtained with the multi-channel 
Reionization And Transients InfraRed camera (RATIR; Butler 
et al. 2012) mounted on the 1 .5 m Harold L. Johnson telescope at 
the Mexican Observatorio Astronomico Nacional on Sierra San 
Pedro Martir in Baja California, Mexico. (Watson et al. 2012). 
Typical observations include a series of 60 s exposures, with 
dithering between exposures. RATIR’s fixed infrared (IR) filters 
cover half of their respective detectors, automatically providing 
off-target IR sky exposures while the target is observed in the 
neighboring filter. Master IR sky frames were created from a 
median stack of off-target images in each IR filter. No off-target 
sky frames were obtained with the optical CCDs, but the small 
galaxy size and sufficient dithering allowed for a sky frame to 
be created from a median stack of all the images in each filter. 
Flat-field frames consist of evening sky exposures. Given the 
lack of a cold shutter in RATIR’s design, IR “dark” frames are 
not available. Laboratory testing, however, confirms that dark 
current is negligible in both IR detectors (Fox et al. 2012). The 
data were reduced, coadded, and analyzed using standard CCD 
and IR processing techniques in IDL and Python, utilizing online 
astrometry programs SExtractor (Bertin & Arnouts 1996) and 
SWarp . 13 Calibration was performed using a single comparison 
star in the SN field that also has reported fluxes in both the Two 
Micron All Sky Survey (Skrutskie et al. 2006) and the SDSS 
Data Release 9 Catalog (Ahn et al. 2012). 

We did not yet possess template images of the host galaxy 
(i.e., prior to the SN discovery or when the SN has faded to 
invisibility) to subtract from the images with the SN present in 


1 1 http://idlastro.gsfc.nasa.gov/contents.html 

12 http://www.sdss.Org/dr7/algorithms/sdssUBVRITiansform.html# 
Lupton2005 

13 SExtractor and SWarp can be accessed from 
http://www.astromatic.net/software. 
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Figure 2. BVRIzJH light curves of SN 2013df from KAIT (solid points) and 
RATIR (open points; with r and i converted to R and /, respectively, following 
http://www.sdss, org/dr7/algorithms/sdssUBVRlTransform,html#Lupton2005). 
The observed curves have been offset from each other for clarity. Shown for 
comparison are the BVRIJH light curves of SN 1993J (Richmond et al. 1994, 
1996; Matthews et al. 2002; dashed lines) and of SN 201 ldh (Arcavi et al. 
2011; Van Dyk et al. 2013b; Ergon et al. 2013; dot-dashed lines), and the 
z-band light curve of SN 2008ax (Pastorello et al. 2008, solid line), all adjusted 
in time and relative brightness. The estimated time t = 0 of the explosion is 
indicated. 

(A color version of this figure is available in the online journal.) 

each band. Consequently, the photometry presented here should 
be considered preliminary; however, the SN is far from the 
main light from the host galaxy, so results including template 
subtraction prior to photometry might not substantially differ 
from what we present here. 

2.2. Early-time Spectroscopy 

We have obtained a number of spectra of the SN at early 
times, using the Kast spectrograph (Miller & Stone 1993) on the 
Lick Observatory 3 m Shane telescope and the DEep Imaging 
Multi-Object Spectrograph (DEIMOS; Faber et al. 2003) on the 
Keck-II 10 m telescope. Ultraviolet (UV) spectroscopy has also 
been obtained using the Hubble Space Telescope ( HST) Space 
Telescope Imaging Spectrograph as part of program GO- 13030 
(PI: A. V. Filippenko). The results of these HST observations 
will be presented in a future paper, together with the bulk of 
the ground-based optical spectra. However, here we present and 
analyze a representative spectrum obtained on 2013 July 11.2 
with DEIMOS using the 600 1 mm~ 1 grating. 

2.3. HST Imaging 

The region of the host galaxy containing the SN site was 
observed with the HST Wide Field Planetary Camera 2 on 1999 
April 29 by program GO-8400 (PI: K. Noll), as part of the 
Hubble Heritage Project. The bands used were F439W (two 
individual images with 40 s exposure times and two with 1 000 s), 
F555W (four 400 s exposures), F606W (two 60 s exposures), 
andF814W (two 40 sand four 400 s exposures). The F555W and 
F814W data were combined with images obtained by programs 
GO-5397 and GO-5972 at an earlier time for the rest of the 



Figure 3. Absolute Flight curve for SN 20 1 3 df( solid points), compared to those 
at VforSNe 1993J (dashed line; colored red in the online version), 2008ax (solid 
line; colored green in the online version), and 201 ldh (dot-dashed line; colored 
blue in the online version) shown in Figure 2. All of the curves are displayed 
relative to the day of V maximum; see the text. 

(A color version of this figure is available in the online journal.) 

galaxy, and drizzled into mosaics in each band at the scale 
ff.'05 pix 1 by Holwerda et al. (2005). 

We have also observed the SN on 20 1 3 July 1 5 with HST using 
the Wide Field Camera 3 (WFC3) UVIS channel in F555W, as 
part of our Target of Opportunity program GO-12888 (PI: S. 
Van Dyk). The observations consisted of 28 exposures of 5 s; 
the short exposure time was intended to avoid saturation in each 
frame by the bright SN. 

3. ANALYSIS 
3.1. Light Curves 

In Figure 2 we display the early-time KAIT and RATIR light 
curves in all bands for SN 2013df. We also include the very 
early V measurement from June 11.202 by Stan Howerton. 14 
We compare these curves with those at BVRIJH for SN 1993J 
(Richmond et al. 1994, 1996; Matthews et al. 2002) and for 
SN 201 ldh (Arcavi et al. 2011; Van Dyk et al. 2013b; Ergon et al. 
2013). We also compare the z-band light curve of SN 2008ax 
(Pastorello et al. 2008). The curves for the comparison SNe were 
adjusted in time and relative brightness to match the curves of 
SN 2013df, particularly at the secondary maximum in each 
band. Clearly, from its overall photometric similarity with the 
other SNe, SN 2013df appears to be an SN lib; see Arcavi et al. 
(2012) for a general description of SN lib light-curve shapes. 
What is most notable from the comparison is that the post-shock- 
breakout cooling of SN 2013df occurred at a later epoch in all 
bands, relative to that of SN 1993J. The post-breakout decline 
of SN 201 ldh occurred at an even earlier relative epoch (e.g., 
Arcavi et al. 2011; Bersten et al. 2012). 

We show in Figure 3 the absolute V-band light curve of 
SN 2013df, relative to those of SNe 1993J, 2008ax, and 201 1 dh. 


14 http://www.flicki-.com/photos/watchingthesky/9035874997 
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(We do not make this comparison based on bolometric or 
pseudo-bolometric luminosity, since, as noted, we consider the 
present photometry of SN 2013df to be preliminary.) The light 
curve of SN 2008ax is from Pastorello et al. (2008), rather than 
from Taubenberger et al. (201 1), since the former has somewhat 
more complete coverage through the secondary maximum. The 
epochs of the secondary maxima in V for SNe 1 993 J, 2008ax, 
and 201 ldh are from Richmond et al. (1996), Pastorello et al. 
(2008), and Ergon et al. (2013), respectively. This epoch for 
SN 2013df was determined by comparing its light curve to 
that of each of the three comparison SNe. The light curves 
have all been adjusted for extinction and for the distances to 
their host galaxies (distance to SN 2013df, above; 3.6 Mpc 
for SN 1993J, Freedman et al. 2001; 8.4 Mpc for SN 201 ldh, 
Van Dyk et al. 2013b; and 9.6 Mpc for SN 2008ax, Pastorello 
et al. 2008). We discuss below our estimate of the extinction 
for SN 2013df. Ergon et al. (2013) have estimated the total 
extinction for SN 1993J and SN 201 ldh to be Ay — 0.53 
and Ay — 0.22 mag, respectively, and Pastorello et al. (2008) 
estimate Ay = 0.93 mag for SN 2008ax (all assuming a Cardelli 
et al. 1989 reddening law and R v — 3.1). 

SN 2013df appears to be the least luminous of all of these 
SNe. The differences in the peak luminosity likely arise from 
differences in the synthesized 56 Ni mass (e.g., Perets et al. 2010). 
For comparison, the 56 Ni produced in the comparison SNe was 
—0.07-0.15 M q , ~0. 10-0.15 M 0 , and -0.06 M 0 in SN 2008ax 
(Pastorello et al. 2008; Taubenberger et al. 2011), SN 1993J 
(Young et al. 1995; Richardson et al. 2006), and SN 201 ldh 
(Bersten et al. 2012), respectively. 

From the overall comparison of the SN 2013df light curves 
in Figures 2 and 3 with the light curves of SNe 1993J and 
201 ldh, we can make an estimate of the explosion date based, 
in particular, on when SN 2013df appears to have reached 
secondary maximum in B VRI. We find that all of the light curves 
indicate that this date is JD 2,456,447.8 ± 0.5, or about June 
4.3 (indicated in Figure 2). This is certainly consistent with the 
earliest discovery epoch of June 7.87 (JD 2,456,451.37). We are 
unable to constrain the explosion date based on KAIT SN search 
monitoring, since the observation of the host galaxy with KAIT 
prior to discovery was on May 25 (i.e., about 10 days before 
discovery), at a limiting (unfiltered) magnitude of 18.5. 

3.2. Spectrum 

In Figure 4 we show the Keck/DEIMOS spectrum of 
SN 2013df and a comparison with spectra of the SNe lib 1993J 
(Filippenko et al. 1993), 2008ax (Pastorello et al. 2008, obtained 
from SUSPECT 15 ; the data are also available at WISeREP 16 ), 
and 201 ldh (available from WISeREP), all at nearly the same 
age. The SN 1993J explosion date was 1993 March 27.5 (Lewis 
et al. 1994), so the spectrum of SN 1993J from April 30 shown 
in the figure is at an age of —34 days. The age of the spectrum of 
SN 2008ax is —30 days (Pastorello et al. 2008). The age of the 
SN 201 ldh spectrum is —28 days (assuming an explosion date 
of 2011 May 31.275; Arcavi et al. 2011). From our estimate 
of the explosion date, above, the SN 2013df spectrum shown 
in Figure 4 was obtained on day —37. This spectrum bears a 
greater similarity with that of SN 1993J than SN 2008ax. Sev- 
eral of the He i lines, notably the 76678 line, atop the broad Ha 
emission profile, appear to have been stronger in the SN 1993J 


15 http://suspect.nhn.ou.edu/~suspect/ 

16 Yaron & Gal- Yam (2012), 

http ://w w w. weizmann .ac.il/ astrophy sics/wiserep/. 
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Figure 4. Optical spectrum of SN 2013df obtained on 2013 July 11.2 
with DE1MOS on the Keck 11 10 m telescope. Also shown for comparison 
are the spectra, at approximately the same age, of SN 1993J from 1993 
April 30 (Filippenko et al. 1993), SN 2008ax from 2008 April 2 (Pastorello 
et al. 2008), and SN 201 ldh from 2011 June 29 (unpublished; from WISeREP, 
Yaron & Gal- Yam 2012). Several spectral features are indicated. 

spectrum (even at a somewhat earlier age) than for SN 2013df. 
This could imply that the H layer for SN 2013df may have been 
more substantial (larger mass, larger radius) at explosion, as 
suggested as well by the post-breakout light-curve behavior of 
SN 2013df, relative to SN 1993J. 

Visible in the spectrum of SN 2013df are absorption fea- 
tures due to Nai D. In Figure 5 we show these features 
after the overall continuum in the spectrum has been normal- 
ized. One pair of lines appears to be weakly visible at effec- 
tively zero redshift, which would correspond to the Galactic 
foreground extinction contribution; these are indicated as “MW” 
in the figure. Another pair is far more distinct, at a red- 
shift of 0.002874, indicated as “Host” in the figure. This is 
consistent with the redshift for NGC 4414, 0.002388, given 
by NED. 17 We have measured the equivalent width (EW) of 
each of the Nai features, D\ and D 2 , for each system to be 
EW(D 2 ) = 0.097 ± 0.011 A and EW (DO = 0.049 ± 0.001 A 
from the Galactic system, for a total EW(Nai) = 0.146 A. 
For the host-galaxy system, EW(D 2 ) = 0.433 ± 0.002 A and 
EW(Dj) = 0.259 ± 0.001 A, for a total EW(Nai) = 0.691 A. 
The ratio EW(D 2 )/EW(Di) is 1.98 for the Galactic component 
and 1.67 for the host-galaxy component. The ratio of the oscil- 
lator strengths of these two lines is 2.0; this intrinsic value is 
approached only at the lowest optical depths (Munari & Zwit- 
ter 1997). We can see that the EW(Nai) [host] is —4.7 times 
that of EW(Nai) [Galactic]. If we assume the Galactic fore- 
ground extinction estimated by Schlafly & Finkbeiner (2011, 
Ay = 0.053 mag), then from the ratio of these two components, 
it is plausible that the extinction internal to the host galaxy is 
Ay = 0.25 mag. 

Poznanski et al. (2011) have stressed that one cannot ac- 
curately infer interstellar visual extinction from EW(Nai) in 


17 The NASA/IPAC Extragalactic Database, http://nedwww.ipac.caltech.edu. 
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Figure 5. Portion of the spectrum of SN 2013df shown in Figure 4, after 
normalization of the continuum, focusing on the region including the interstellar 
Nai D AA.5890.0, 5895.9 absorption features. The possible feature associated 
with the Galactic component is indicated by the dashed lines as “MW,” and 
the feature (corresponding to a redshift of 0.002874) associated with the host 
galaxy is indicated as “Host.” 

low-resolution spectra. The resolution of the Keck DEIMOS 
spectrum may be just at the limit, however, where we can use 
the more accurate relations from Poznanski et al. (2012). If 
we apply their Equation (7) for EWl/T), we obtain EiB — 
V) = 0.11 ± 0.04 mag; from their Equation (8) for EW(D0, 
we obtain E{B — V) = 0.08 ± 0.04 mag; and from their 
Equation (9) for the total EW(£>i + Di), E(B — V) = 0.09 ± 
0.01 mag. The uncertainty-weighted mean from all three esti- 
mates is then E(B — V) — 0.09 ± 0.01 mag. If we assume the 
Cardelli et al. (1989) reddening law, with R v = 3.1, then the 
host-galaxy extinction is A v = 0.28 ± 0.04 mag. (We note that 
Poznanski et al. 2012 assume the older, somewhat higher values 
of Galactic extinction from Schlegel et al. 1998.) This estimate 
of the host-galaxy extinction is quite consistent with the one we 
made above, which we adopt. We therefore assume hereafter 
that the total (Galactic foreground plus host galaxy) extinction 
for SN 2013df is Ay = 0.30 mag. We adopt an uncertainty in 
the extinction of ±0.05 mag. 

3.3. Progenitor 

We display in Figures 6(a) and 6(b) the subregion of the 
WFPC2 mosaic around the SN position at F555W and F814W, 
respectively. We show the combination of all the WFC3 expo- 
sures in Figure 6(c). From a precise comparison between the 
pre-explosion WFPC2 images and post-explosion WFC3 im- 
ages, we can identify the SN progenitor. We had previously 
attempted an identification of the star (Van Dyk et al. 2013a), 
but this was employing much lower-resolution, ground-based 
images of the SN. None of the three objects that we had pre- 
viously nominated as candidates, numbered in Figure 6, turned 
out to be the likely SN progenitor. Using 1 1 stars in common 
between the WFPC2 and WFC3 images, we astrometrically reg- 
istered the images to an accuracy of 0. 1 1 drizzled WFPC2 pixel, 


Table 1 

HST Photometry of the SN 2013df Progenitor 3 


F439W b 

fi b 

F555W 

V 

F814W 

/ 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

>25.65 

>25.54 

24.535(071) 

24.515 

23.144(055) 

23.106 


Notes. 

a Uncertainties (lcr) are given in parentheses as millimagnitudes. 
b 3tr upper limit. 


or 5.5 mas, and found, instead, that the SN corresponds to the 
position of the star indicated by the tick marks in Figures 6(a) 
and 6(b). We therefore identify this star as the likely progenitor 
of SN 2013df. 

We extracted photometry for the source from the pre- 
explosion WFPC2 images using Dolphot v2.0 (Dolphin 2000). 
(We disregarded the short F439W, F606W, and F814W ex- 
posures, since the signal-to-noise ratio is quite low in these.) 
The progenitor location is in the WFPC2 chip 3. The Dolphot 
output indicates that the progenitor is likely stellar, with an 
“object type” of 1; the “sharpness” parameter, at —0.37, is 
slightly outside the acceptable range (—0.3 to +0.3) for a “good 
star” (Dolphin 2000; Leonard et al. 2008). However, if we run 
Dolphot on the F555W and F814W images only (excluding the 
F439W images), this parameter has a more acceptable value 
of —0.25. We therefore consider it most likely that this is a 
well-resolved stellar object. The star is not detected by Dolphot 
in the F439W images (nor is it detectable upon visual inspec- 
tion of the images). However, we estimate a 3er upper limit 
to its detection. The Dolphot output automatically includes the 
transformation from HST flight-system magnitudes to the corre- 
sponding Johnson-Cousins magnitudes (Holtzman et al. 1995). 
We present the Dolphot results for the progenitor in Table 1 . We 
also show the resulting spectral energy distribution (SED) of the 
star in Figure 7. 

To our knowledge, no measurement exists of the metallicity 
in the vicinity of the SN position (such as from spectroscopy of 
nearby Hu regions). So, we are unable to quantify accurately 
the metallicity at the SN site with available data. However, 
we can provide at least an approximate estimate, based on the 
assumption that the abundance gradient in the host spiral galaxy 
follows the same behavior as for the sample of spirals analyzed 
by Zaritsky et al. (1994). After deprojecting NGC 4414 from 
its inclination and position angle (56° and 160°, respectively; 
Jarrett et al. 2000), we find that SN 2013df occurred ~31" (or 
~2.5 kpc) from the nucleus. From the Zaritsky et al. nuclear 
abundance and gradient, then, the oxygen abundance (a proxy 
for metallicity) at the SN site is 12 + log [O/H] ~ 8.8, which is 
comparable to the solar abundance (8.69 ± 0.05; Asplund et al. 
2009). Therefore, it is reasonable to assume that the SN site is 
of roughly solar metallicity. 

We compared the observed SED of the progenitor with 
synthetic SEDs derived from MARCS (Gustafsson et al. 2008) 
supergiant model stellar atmospheres. 18 The model atmospheres 
are spherical, with standard composition at solar metallicity, 
surface gravity log g = 1.0, and microturbulence of 5 km s _1 . 
The spherical-geometry model atmospheres were computed, 
generally, for stars with masses 0.5, 1.0, 2.0, 5.0, and 15 M 0 ; 
we chose models for this last mass as most appropriate for 
the likely massive progenitor, given what has been inferred 
for the progenitor initial masses for other SNe lib. We found 

18 http://marcs.astro.uu.se/ 
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Figure 6. Portion of the archival HST WFPC2 images of NGC 4414 from 1999 in (a) F555W and (b) F814W. The likely progenitor of SN 2013df is indicated by tick 
marks, (c) A portion of the HST WFC3 F555W image of SN 201 3df, to the same scale and orientation. The SN is indicated by tick marks. We have astrometrically 
registered the WFPC2 and WFC3 images to an accuracy of 5.5 mas. The three progenitor candidates initially identified by Van Dyk et al. (2013a) are indicated (“1,” 
“2,” and "3”) in panel (a). North is up. east is to the left. 


that the synthetic V — I colors , in particular, derived from models 
at a given temperature are essentially mass-independent (they 
differ by ~0.03 mag). The models were reddened based on 
the assumed total extinction for SN 2013df, above, assuming 
the Cardelli et al. (1989) reddening law, and were normalized 
at V. The model that provided the best comparison with the 
photometry has 7' c ff = 4250 K, which we show in Figure 7. 
Warmer and cooler models did not compare well with the 
observations. We therefore adopt this effective temperature for 
the progenitor, with a generous uncertainty of ±100 K. 

The absolute V magnitude of the progenitor, corrected for the 
assumed distance modulus and extinction to the SN, is My = 
—6.89 ± 0.10. The uncertainty arises from the uncertainties in 
the photometry, the extinction, and the distance modulus, all 
added in quadrature. The uncertainty in the transformation from 
flight system to Johnson-Cousins magnitudes is <±0.01 mag 
(Holtzman et al. 1 995). The bolometric correction obtained from 
the MARCS model at 7’ c n = 4250 K is BCy — —0.71 mag. 
The uncertainty in BCy, resulting from the uncertainty in 


T e ff, is ~0.05 mag. Including this uncertainty, the bolometric 
magnitude is Mb 0 i = —7.60 ± 0.15. Assuming a bolometric 
magnitude for the Sun of 4.74, this corresponds to a luminosity 
L bol = 10 4 - 94±0 06 L 0 . We show the locus for the progenitor in 
a Hertzsprung-Russell (H-R) diagram in Figure 8. 

4. DISCUSSION AND CONCLUSIONS 

The SN 2013df progenitor’s position in the H-R diagram is 
significantly blueward of the RSG terminus of model massive- 
star evolutionary tracks with rotation, such as the 15 Mq track 
(Ekstrom et al. 2012) in Figure 8. This indicates that the 
star likely has had its envelope somewhat stripped by some 
mechanism. The expectation so far for the progenitor scenario 
for SNe lib is an interacting massive star binary system, such 
as for SN 1993J (e.g., Podsiadlowski et al. 1993; Maund et al. 
2004) and SN 201 ldh (Bersten et al. 2012; Benvenuto et al. 
2013). (However, see Crockett et al. 2008 for other possibilities 
regarding the progenitor of SN 2008ax.) A binary channel is also 
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Figure 7. Spectral energy distribution of the probable progenitor star of 
SN 20 1 3df, based on the photometry presented in Table 1 , Shown for comparison 
is synthetic photometry derived from a MARCS (Gustafsson et al. 2008) 
supergiant model atmosphere at solar metallicity with effective temperature 
4250 K, reddened assuming Ay = 0.30 mag and the Cardelli et al. (1989) 
reddening law. The model has been normalized at V. The hashed region in 
the diagram represents the range in brightness for the model within the total 
uncertainty in the observed V magnitude. 
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Figure 8. Hertzsprung-Russell diagram, showing the locus of the probable 
progenitor of SN 2013df (solid symbol; colored red in the online version). For 
comparison we also illustrate the massive-star evolutionary tracks with rotation 
from Ekstrom et al. (2012) at initial masses 12, 15, and 20 Mq (curves), as well 
as the locus of the SN 1993J progenitor (dashed symbol; colored blue in the 
online version) from Maund et al. (2004). 

(A color version of this figure is available in the online journal.) 


strongly favored for SNe lib, in general, based on statistical 
arguments accounting for their observed relative rate among 
core-collapse SNe (Smith et al. 2011), as well as on the small 
ejecta masses and H/He line ratios compared to detailed models 
(Dessart et al. 2011; Hachinger et al. 2012). A blue binary 
companion is expected to survive the explosion for SNe lib 
(Maund et al. 2004; Benvenuto et al. 2013). The locus of the 
SN 2013df progenitor in the H-R diagram is quite similar to, 
albeit possibly somewhat less luminous than (though within the 
uncertainties), that of the supergiant progenitor of SN 1993J 
(Maund et al. 2004, see also Figure 8). 

We can attempt to assign an initial mass to the progenitor, 
which, following Van Dyk et al. (2011) for SN 201 ldh, would 
have been ~16 — 17 M 0 , based on comparing the locus in T e g and 
Lboi to the corresponding values of the 15 and 20 M 0 tracks in 
the figure. On the other hand, adopting the approach of Maund 
et al. (2011) of basing the likely initial mass on the luminosity 
at the terminus of the available tracks, we estimate the mass 
to be ~13— 14 M q . As Maund et al. (and references therein) 
pointed out, the final luminosity at the top of the RSG branch 
is more relevant to the condition of the yellow supergiant. The 
timescale for the transition redward for massive stars across 
the H-R diagram is comparatively fast, and the stars will 
subsequently sit at the RSG locus for the majority of their post- 
main-sequence evolution. In a binary model, mass transfer will 
truncate the outer radius of the donor star, resulting in a hotter 
photosphere; however, the core nuclear-burning luminosity will 
be relatively unaffected. Note that, if the progenitor experienced 
mass transfer in a binary, substantial mass stripping can lower 
the core temperature and, hence, lower the surface luminosity. 
So, the progenitor could have been somewhat more massive 
initially, perhaps ~15 Mq. 


The effective radius of the SN 2013df progenitor 14 yr 
before explosion, based on our estimates of T e g and Lb 0 i, was 
/( c it = 545 ± 65 Rq. Unfortunately, we cannot produce an 
independent estimate from the light curves using the relations 
between SN progenitor radius and photospheric temperature 
by Chevalier & Fransson (2008), Nakar & Sari (2010), and 
Rabinak & Waxman (2011), since these relations tend to break 
down only a few days after explosion, and the earliest set of 
photometric data points for SN 201 3df is at an age of ~10 days. 
This star was evidently more extended than the progenitor of 
SN 201 ldh — the yellow supergiant progenitor of SN 201 ldh 
modeled by Bersten et al. (2012) has radius R ~ 200 R e . The 
later post-shock breakout cooling for SN 2013df implies that its 
progenitor may have had a larger radius than that of the SN 1 993 J 
progenitor. However, both Maund et al. (2004) and Van Dyk 
et al. (2013b) calculated that the progenitor of SN 1993J had a 
radius of ~600 R a , which is somewhat greater than, although 
within the uncertainties of, what we have estimated for the 
SN 2013df progenitor. In order for the SN 1993J progenitor to 
have had a smaller effective radius, it must have been hotter 
(>4265 K) and of lower luminosity (<10 5 1 L a ) than currently 
understood (e.g., Maund et al. 2004). To have a larger radius, 
the SN 2013df progenitor would have had to be cooler and also 
potentially more luminous. However, to be cooler would require 
the extinction to be even lower than the small amount that we 
have estimated. Additionally, the somewhat larger bolometric 
correction with the lower temperature would be offset by the 
smaller extinction correction, leading, ultimately, to very little 
change in the estimated luminosity. At the least, we can say that 
the SN 2013df and SN 1993J progenitor stars appear to have 
been relatively comparable in nature. 

In summary, we have shown that SN 2013df in NGC 4414 
is an SN lib, based on early-time photometry and 
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spectroscopy, and that its properties are most similar to those of 
SN 1 993J. However, the mass of 56 Ni synthesized in the explo- 
sion likely differed from that in SNe 1993J, 2008ax, and 201 ldh. 
Furthermore, the light curves of SN 2013df — specifically, the 
late cooling from the post-shock-breakout maximum — indicate 
that its progenitor must have been more extended in radius than 
that of SN 1993J. Using archival pre-SN, high spatial resolution 
HST images, together with HST images of SN 201 3df, we have 
identified a star that we consider to be the likely SN progen- 
itor; it had T eB ~ 4250 K and Lb 0 i ~ 10 4 ' 94 L 0 . Ultimately, 
verification of this identification of the SN 2013df progenitor 
should occur when the SN itself has greatly faded, as has been 
done for the SNe lib 1993J (Maund & Smartt 2009) and 20 1 ldh 
(Ergon et al. 2013; Van Dyk et al. 2013b). This can be accom- 
plished with high-quality images obtained with HST. Deep, very 
late-time HST imaging in the blue or ultraviolet could reveal a 
putative blue binary companion star at the SN location, although 
the fact that SN 2013df is at least a factor of two more distant 
than both SNe 201 ldh and 1993J could make such a detection 
a more challenging prospect. 
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